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We generalize the consistency condition for the three-point function in single field 
inflation to the case of dissipative, multi-field, single-clock models. We use the re- 
cently introduced extension of the effective field theory of inflation that accounts 
for dissipative effects, to provide an explicit proof to leading (non-trivial) order in 
the generalized slow roll parameters and mixing with gravity scales. Our results 
illustrate the conditions necessary for the validity of the consistency relation in situ- 
ations with many degrees of freedom relevant during inflation, namely that there is 
a preferred clock. Departures from this condition in forthcoming experiments would 
rule out not only single field but also a large class of multi-field models. 
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1. INTRODUCTION 

The squeezed limit of the three-point function of curvature perturbations, geometrically: 

((k L (k s (k s )k L ^0 <-> ^ f kh — > , 

k s 

is tailor made to probe the dynamics of the early universe, due to its intrinsic dependence on 
the very mechanism that produces such perturbations. It is observationally relevant in many 
contexts, for example in measurements of large scale structure or CMB /i-distortion 

(3)]. For models with only one light dynamical field that is relevant during inflation, there 



is a consistency condition which relates the three-point 
deviations from scale invariance in the power spectrum 
more precisely as 



unction in the squeezed limit to 



Mi. 



This relation can be written 



hm (C kl CM = -(27r) 3 5 3 ^ P c (k L )P c (k s ) [(n s - 1) + O (kl/k 2 s )} , (1) 

where k x — > k^ and (k 2 — k 3 )/2 — > ks, for the long- and short-wavelength modes respectively. 
Also P({ki) is the power spectrum, i.e. (CkXkj) — ^ 3 (ki + kj)P^(ki), and n s — 1 is the tilt 

_ rflogfc 3 P c (fc) 
Hs L ~ dlogk ' {Z) 

such that (CC) — k~ z+ris ~ l . Alternatively, defining F(ki, k 2 , k^) by 

(C kl CM = (2vr) 3 5 3 (^k l )F(fc 1 ,A; 2 ,A ; 3), (3) 

i 

and using the definition of the parameter /^l 

fSq = 5 F(k 1 ,k 2 ,k 3 ) 

" *s-o 6 (P c (h)P c (k 2 ) + Pdh)P C (k3) + P c (k 2 )P c (h)y 1 ] 

the consistency relation in the exact squeezed limit can be written as 

/n! = ~{n. - 1). (5) 
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The expression in Eq. ([T]) follows in the long- wavelength limit, k^ — > 0, as an expansion 
in powers of ki/ks which is valid to all orders in the slow roll parameters js-8|. Therefore, 
Eq. ([I]) would remain valid regardless of the value for n s — 1, even for hypothetical scenarios 
with large departures from scale invariance. 

Eq. ([I]) is nothing but the manifestation of the fact that the long mode has no dynamical 
effect on the short modes, and that the three-point function in the squeezed limit is in a 
sense just an artifact of the coordinates being used. If rather than co-moving coordinates 
we perform a transformation to a more physical set of coordinates, i.e. x l phys = e^ L x l , or 
fcp hys = e~^ L ks then the amplitude of the small scale power would become independent of 
the long mode. There is no three-point function in these coordinates. In gravity physical 
effects are encoded in the curvature and thus one expects that the first correction in Eq. (pQ) 
scales as (ki/ks) 2 - More detailed versions of this argument can be found in 5|48]. 

The validity of the consistency relation is evident in cases where heavy (or decoupled) 
fields are present, but becomes non-trivial in situations where light excitations are pro- 
duced during inflation which subsequently decouple, or in general when additional degrees 
of freedom (ADOF) produce negligible direct contributions to (, but significantly alter the 
dynamics of the (single) clock which controls the end of inflation. 

Scenarios where ADOF couple to the clock, but do not contribute to density perturba- 
tions, are tantamount to study dissipative effects during inflation, such as warm [9] and 
Tapped 10| inflation, which were recently incorporated in 11] (based on ideas originated in 
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13]) within the effective field theory (EFT) of inflation formalism developed in 14J-I23]. 



An explicit proof of the consistency condition for single field inflation, valid at leading 
non-trivial order in the slow roll parameters, was given in [6| using the EFT approach 
[yj]. Our purpose in this paper is to extend the results in js] to single-clock models with 
dissipative ADOF. Proving the consistency condition in these type of scenarios is more 
challenging than in single field inflation due to the presence of ADOF whose detailed 
dynamics, and couplings, remain vastly unknown. However, thanks to the EFT approach 
fill . Q], we will show that under certain circumstances the consistency relation holds to 
leading (non-trivial) order in the generalized slow roll parameters and mixing with gravity 
scales. Departures from this condition in forthcoming experiments would rule out not only 
single field, but also a large class of multi-field single-clock models. 
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This paper is organized as follows. In the next section we review the EFT formalism 



introduced in 
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14| , in particular: The assumptions we make about the properties of the 
noise and response for the ADOF; the constraint equations in the presence of ADOF and 
the constancy of ( outside the horizon; the computation of the power spectrum and tilt. 
Then, in sec. |3]we explicitly show the validity of the consistency relation to leading (non- 
trivial) order in the generalized slow roll parameters, first in the limit M p — > oo and later on 
including the mixing with gravity. As expected, in sec. H]we show that even in the presence 
of ADOF a long-wavelength mode amounts to a rescaling of the coordinates, as in single 
field inflation. 

We concentrate in models with scalar ADOF, although we discuss towards the end in sec. 
|5] how to generalize our analysis for vector and tensor operators. We relegate details and 
more technical points to appendices, including the different contributions from each source 
term and {ki/ks) 2 scaling in the three-point function. Everywhere we set c = H = 1 and 
adopt the mostly plus sign convention. 

2. EFFECTIVE FIELD THEORY SET UP 



Here we introduce the necessary elements of the EFT formalism developed in 11( to study 
dissipative effects during inflation. 1 We concentrate on the features which are relevant to 
analyze the squeezed limit of the three-point function. 



2.1. The effective action 



In UJ dissipative effects are incorporated via the coupling to a set of (scalar, vector, 
tensor) composite operators Q^ v - in the effective action, which is constrained solely by 
invariance under all the symmetries of the long distance physics. 

The so-called unitary gauge is defined such that the perturbations of the preferred clock 
vanish. In this gauge the unit vector perpendicular to surfaces of constant time t takes 
the form n p = — 5°(— g 00 )^ 1 ^ 2 , and the extrinsic curvature of the surfaces is = g w V p n v , 
where g w = g pp +n p n p is the induced spatial metric. In addition, one can choose coordinates 



1 The EFT of multi- field inflation with many light degrees of freedom was studied in [24[ . 
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so that the metric for fluctuations around the quasi-deSitter background is given by 

ds 2 = -N 2 dt 2 + a 2 {t)5 ij e 2<: {dx i + N i dt)(dx j + N j dt). (6) 
In this paper we ignore tensor perturbations. 

The action in the unitary gauge reads 

M 2 r i r 

S = -f J d*Xyf^R+-J d 4 x^(p-p-(p + p)g 00 ) 

+ \J d'x^M^m + g 00 ) 2 -^ I d 4 x^Ml(t)SKji(l + g 00 ) 

- [ d 4 x^f(t)0 + S + --- , (7) 



where M 2 = (StiGn) -1 , and O is a scalar composite operator associated with the dissipative 
degrees of freedom and So represents the action for this sector, which we do not need to 
specify. The dots stand for higher derivative corrections. Also, bared quantities denote 
their unperturbed value. Throughout the paper we assume the presence of an approximate 
shift symmetry, such that functions of time appearing in the action change very little in a 
Hubble time. This will be the basis of our generalized slow roll approximation. Since we 
are only interested in the squeezed limit of the three-point function the higher derivative 
terms are irrelevant. 



The field ir, which characterizes the perturbations of the preferred clock, is then intro- 
duced via the Stiickelberg's trick, i.e. t — y t + tt [14]. It is easy to show that the relation 
between 7r and ( is the same as in single field inflation. Taking into account that in the 
7r-gauge the metric takes the form 



ds 2 = -N 2 dt 2 + a 2 (t)8 ij (dx i + N l dt){dx j + N j dt), 



(8) 
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the relationship becomes, up to second order in the perturbations [4|, 

C = -Hit + Htttt + ^ifvr 2 + a, (9) 
a = ^2 {-di-Kdiix + d^didjidiirdjn)) . (10) 

Notice we restricted ourselves to scalar couplings for the ADOF of the form f(t)0, which 
naturally produce a 777 dissipative term when the Green's function obeys ImG (cu) ~ u; 
fill ]. There are, however, other type of scalar couplings, such as Og 00 . Even though these 
produce terms involving derivatives of tt, they can also contribute at leading order in the 
slow parameters in cases when ImG (u) ~ as it is required for 77V dissipation 11] . We 



will return to this possibility, as well as having vector and tensor couplings, later on in sec. |5j 



2.2. Noise &: Response 

Notice that, neglecting metric perturbations, the presence of the ADOF affects the clock 
only through the coupling between tt and O. As explained in [ll], this introduces dissipation 
and noise effects into the dynamics of tt which can be characterized by the splitting of O 
into pieces: 2 

0(t, x) = 0{t) + 60 s (t, x) + dO R (t, x), (11) 

where 0(t) is the background value, 50fi(t,x) represents the change in the expectation 
value that results as a response to the interaction with tt, and SOs{t,x) corresponds to a 
stochastic noise. We will work in the approximation where both the response part 50n(t, x) 
and the two-point correlator of the noise (SOs(t, x)80s(t', x')) are local in time (namely we 
neglect any memory effect), and we assume an expansion in spatial derivatives holds. Then, 
in unitary gauge, we have for the noise 

(80s(t,k)SO s (t',q)) ~ ^^-t')(2vr) 3 5 (3) (q + k). (12) 



2 We use SO s , instead of 80s, m order to avoid confusion with 5 being used to describe short modes. 
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On the other hand, for the response part we expect in general a (local) relationship of the 
sort [ll | 

O + SOI ~ A (t) F [y/=j?»,K*,t] , (13) 
where F is some generic function and Aq provides the necessary scale such that 



6(t) = A (t)F[l,3H,t}. 



(14) 



Because of the time dependence of the background, it is convenient to define the operator 
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SOftiir) = (O r - 6{t + tt)) 



(15) 



so that going from the unitary to the 7r-gauge we have SO R — > 50 s r (tt). (G is a scalar 
operator which does not introduce any 7r's. This replacement is then required to cancel out 
the expansion in t + tt from the background value 0(t).) 
Using the expression in Eq. ffl3l) we have 



50 s R {u) = A a (t + ix)F (rfd^t + tt), K£(t + tt), t + tt) - 0(t + tt) 



(16) 



where 



-g^d u (t + n) 



(17) 



y/-g v Pd v (t + 7r)d p (t + n) 

Using analyticity for small perturbation around the background n M = 5^°, Eq. (TIB"]) can 
be expanded as 3 



-V (t + tt) 



'l+g^d,(t + 7i)d u (t + n)) + 



(18) 



Wait + tt)(1 + gTdp{t + 7r)9„(t + tt)) 2 + •••]. 



The ellipses in Eq. ( TT8l) account also for terms proportional to the extrinsic curvature 
which one can show at most contribute to a renormalization of the speed of sound, and 
therefore does not alter the proof of the consistency condition. This expression encodes 



For example let us consider a function Or — F(Ai\/X), where X = —g liU d li (t + ir)d l ,(t + , K). Then, X = 1 
and F(A X (t + tt) ^1 + SX) ~ 6{t + tt) + Aj (i + tt)5X + • • • , with SX = -(l+g^d^t + ^d^t + n)). 

Hence we obtain Eq. (jT5J) for the expansion of 50 S R . 
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all the information about the coupling between O and 7r, including all possible non-linear 
interactions. 4 Notice that the RHS of Eq. f lT8|) is (approximately) invariant under tt — > n + c 
provided Vo{t), Wo(t), . . . are slowly varying functions of time. From now on we will assume 



that is the case for all time dependent functions in the problem, including 0(t). 



what we refer as the existence of an approximate (emergent) shift symmetry in tc 11]. 



'his is 



Note also that we write the metric components in the contractions explicitly. This is 
because, as we will discuss momentarily, they play a crucial role to guarantee the constancy 
of ( outside the horizon, as well as to incorporate the mixing with gravity. 



In what follows we will assume an expansion as in Eq. (fT8|) applies, not only to O, but 
all operators that derive from the ADOF, e.g. (Tq)o etc. 



2.3. Constraint equations & Constancy of ( outside the horizon 

Using the expression in Eq. ([7]) the background equations are given by: 

3H 2 M*=p + p o + f(t)0, (19) 
p + 3H(p + p) + fO = } (20) 
Po + 3H(p o +p o ) + fd = 0, (21) 

where the last one is obtained by combining the second and the time derivative of the first, 
and we used that the background stress tensor for the ADOF obtained from So takes the 
perfect fluid form 

T%, = (po + Po)n^n v + g^po- (22) 



4 In other words, Eq. (fT8|) includes both: non-linear response applied to linear forces, as well as linear 
response applied to non-linear forces in 7r. 
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To linear order in the perturbations the constraint equations become 
Ml 



2HMpdi5N + -f {didjNi - dANj) - (p + p)d^ - M\di{8N - tt) + tfsf 5 = 0, (23) 
—6H 2 5NMp - 2Hd i N i M 2 + ^ iV* + 3#(5iV - ei/vr) + ^ J 

+{N C + 3HMl)(5N --k)-pK- fOir - f50 {1) + 6 A® = 0, (24) 

where the superscript ^ indicates the order in the perturbations, and 5A^ and SB^ are 
the linearized part of 

A = jaw' < 25 > 

Note that for the background we have A = —po and = 0. We also introduced 

iV c ^(p + p + 4M 2 4 ), (27) 
which will play the role of our normalization scale for tt later on. 

Solving for the lapse and shift we find (for H 2 M 2 > N c , HM\) 

^ll^-jgj-^-^), (28) 

2HM% 2EM 2 V 2HMI ' 1 ; 

diN* ~ -— ^ (W - etfvr) - 3#5£« _ + - SA^ , (29) 
where 6B<V = d^drfB^ and we used Eq. (EDI. 



It is now convenient to rewrite these equations in terms of 50 s , as in Eq. f JT5|) . and 
accordingly: 

cL4 s = (5A + p vr, (30) 

5B S = 5B + (Po+Po)k- (31) 
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The reason will be soon clear, but notice that at linear order these operators correspond to 
the (scalar) perturbations of O, (Tq) q and d~ 2 di {T®) , respectively. 5 Using these defini- 
tions, and Eq. ( )2T]) . we can then write 

M.^A-^-f), (32) 
2HM* 2HM 2 p ' v ; 

diNi - ~2ikp { Nc{ir " eHll) ~ 3H5B u " ^iJr + f6 ° s ^ ~ M «) • (33) 

The expression in Eq. ( )32|) closely resembles the relationship in single field inflation, 
except for the appearance of new terms proportional to perturbations of the ADOF. However, 
as we argue in the previous section, the response for this term must obey a local expansion 
similar to the one in Eq. ffTS]) . 6 In other words, including the mixing with gravity to leading 
(non-trivial) order, 

8B s {1) ~V B (t)(>k-5N) + --- , (34) 



such that (M% = HV B ) 



Hence, 



fi^ tg , + -»^.(*-a y )-Mfr-f">. (35) 



(5N-eHn) ~ -^—(n-eHn) + -- - , (36) 
1 + £b 



where we defined e B = ^i H i , % = 2 h*m* > so tnat e = e B -e Si- 



V V 



The expression in Eq. ( 1361) now mimics what occurs in single field inflation. In that case 
the second term can be shown to vanish for long- wavelength modes since £l ~ — > ttl ~ 
eHu^. In order to use a similar reasoning we need to show that -^(Hir) — > in the presence 
of the ADOF. Or, in other words, that ( is still conserved outside the horizon. But using 



5 The reader may be more familiar with gauge transformations of the fluid variables associated to the 
linearized stress tensor of a fluid, for which Tq = —5p, d~ 2 d{T® = (p + p)Su. 

6 Here we ignore the extrinsic curvature contribution since 6K£ ~ ec^C + 0(£ 2 ), and therefore it does not 
alter our results to the order wc work here. 
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Eq. (|18[) we have (including the mixing term) 

SO s R ~ Vo(t)(n - 8N) + ■ ■ ■ ~ Vb(t) (l + J (vr - eHn) + ■■■ , (37) 

which is proportional to (, thus it turns off when ( is constant. This implies a constant 
£ = — i?7r remains a solution. Therefore 

tt l ~ eHn L (38) 

in the long- wavelength limit (see appendix |A] for more details). Moreover, from here we 
conclude 

6N L ~ e#vr L . (39) 
A similar argument can be used for Eq. fl33|) . 



Note that in principle not only the response, but also the noise will enter in the above 
expressions for the lapse and shift. However (513$) = 0, so that the noise part only affects 
7r, and C, as a stochastic source and it will not alter its constancy outside the horizon. This 
is because, within the local approximation and for slowly varying functions of time, the 
computation of the power spectrum depends on integrals of Green's functions which are 
solutions of the homogeneous equation away from the 5-source. As long as this equation 
admits ( = constant as a solution, under our assumptions the noise will not modify such 
behavior. 



Notice that in the previous analysis we did not require e to be small, then the constancy 
of ( also applies to FRW cosmologies. However, since we assume the ADOF do not modify 
significantly the background dynamics, we do require 6b <C 1. More generally, by assump- 
tion, the response parts for the ADOF are functions of the geometry, namely g 00 , Kfi, and 
these tend to their unperturbed value when e# — ?■ 0. Therefore -in the case of stable per- 
tubations in the absence of ADOF- this property will not be modified by adding a small 
correction to the background equations. 

One may wonder about the case when 6b — 1. In that scenario the coupling to gravity 
in the (9-sector cannot be treated perturbatively. However, following similar steps as above, 
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one can still show that a constant ( will stay a solution, although we cannot guarantee it will 
be the dominant one 22j . Note this must be the case because £ — > £ + A is a symmetry in an 



FRW background at linear order, since A can be absorbed into a re-scaling of the coordinates. 
Moreover, we expect this to remain valid under the assumption of local dynamics; since the 
response SB R will indeed turn off when ( = in the ki — > limit, and there will be no 
contribution to the long-wavelength curvature perturbations from earlier times when the 
modes were inside the horizon. 

One could still ask if — 1 is even feasible during inflation, namely having an expanding 
background together with a scale invariant spectrum with constant curvature perturbations 
outside the horizon. 7 First of all, the slow roll condition (-H/H 2 ) <C 1 might be violated 
by contributions from the energy density in the (9-sector. We could conceivable devise a 
scenario where violations of slow roll occur during a short period of time. But that is not 
the spirit of most dissipative models, where the energy scales associated with the ADOF are 
smaller than M 2 H. We will not discuss modifications of this set up in this paper, and from 
now on assume es <C 1. 



2.4. Power spectrum & tilt 

We now move on to the computation of the power spectrum, putting emphasis on all 
possible contributions to the tilt. Using Eq. (j7|) and expanding to quadratic order, the 
action for ir can be written as (ignoring higher derivative corrections) 

S* = J d 4 xa 3 ^ jvr 2 - - i-vr (80 R + 50 s ) + • • • j , (40) 

where 



and N c was introduced in Eq 



2 (p + p + HM, 



3 s 

(41) 



s (P + P + 4M!)' 
( )27l) . which can be also written as N c = (p + ~p + HM^/c 2 . 



As we mentioned, we assume the presence of an approximate shift symmetry such 



7 Note that in this case, however, the explicit proof of the consistency condition becomes more involved, 
since one cannot ignore the contribution from the ADOF to SN for short wavelength modes (see sec. I3.2j) . 
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that functions of time appearing in the action change very little in a Hubble time (e.g. 
^ < 1, j^ff <C 1, etc), and also ignored terms proportional to the generalized slow roll 
parameters which do not alter the computation of the power spectrum to the order we 
work here. (We can then replace 50 s R (ir) by 50 R (ii), since O is suppressed by a slow roll 
parameter. We omit the superscripts below.) 



Recall also that, due to the presence of ADOF, our field ir differs from the choice made 
in [14|, denoted as n in [11]. The quadratic Lagrangian for n is uniquely determined 
by the quasi-deSitter geometry, e.g. H, H, whereas in Eq. fT40l) the normalization scale 



N c is a free parameter. (One can show c 2 s N c < —2M^H - assuming the contribution 
from Mi is subleading - if the stress energy tensor of the ADOF obeys the null en- 
ergy condition [11].) Our choice of unitary gauge, however, guarantees the relationship 
between ( and 7r remains given by Eq. Note that this relation is altered if writ- 

ten in terms of 7r since (schematically) 7T ~ n + SO. See the discussion in ll| for more details. 



The equation for tt that follows from Eq. f|40l) becomes 

/ 2 2 -f i 

7t k + 3H7T k + -^n k + -L60$ {t, k) = -jj-SO s (t, k), (42) 

where SO^ is the response part at leading order in tt which, using Eq. ( JT8|) . can be written 

as 

SO^ ~ V {t)rc, with V (t)=jN c /f, (43) 

such that 

/ 2 2 £ 

n k + (3H + 7 )7T fc + -^tt, = -{-50 s (t, k). (44) 

The expression in Eq. (H3l) makes manifest the presence of an approximate (emergent) 
shift symmetry at the level of the response so that, to leading order, we can neglect terms 
which do not involve derivatives. To compute the power spectrum we also need the two-point 
function for the noise which is given in Eq. (fl2]) . 

The solution of Eq. (j44D to linear order is given by (neglecting the homogenous solution 



which dies off in the 7 ^> H limit) [11 1 



711(77, k) 



fkc s 



dvj 'G~{kc 8 \r}\, kc s \r]' 



/in 5 °£ 



(kc s r}') 



where 7/0 is an early enough initial time, and 



G 7 (z,z') = \z [Y v {z)J v {z') - J v {z)Y v {z')\ 



with v = I + 2jEx, 2; = —kc s r), and 2/ = —kc s r)'. Then, from Eq. (jl2jl we obtain 



P 7r (fc) = (7r 1 (r7,k)7r 1 (r / ,k)) 



iV 6 2 (fcc s ) 3 



For fcc s 77 — )■ and kc s rj — »■ —00, and using Pa = H 2 P n , we find 



A ? = A; 3 P ? (fc) 



2„*3 



7rr(^ + 4) 



where the * means that the quantity is evaluated at freeze out c s fc/a(i*) ~ \f%^FT i< 
particular, in the strong dissipative regime, where 7 3> H, we obtain 



H 2 



2c*(c*N*Y 



Therefore, the tilt n s — 1 is given by 



— 1 



rflogA f 1 d 
dloe-A; ~~ IFdF 



log 



^/M^(- + i)3r(^) 4 \ 



iVfcf 



vrr(^ + 4) 



+ 2e} - 3s* - 2(e* Nc + e*) + R ( 7 */#*) (e* + e 



7/) 



where we introduced the generalized slow roll parameters: 

H i>p f 7 c s N c 

€ ~ # 2 ' €uo ~ Hue' C/ ~ Hf' 67 ~ 2HY S ~ Hc s ' " Nc ~ HN C ' 
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and denned the function 



R{x) = \x \-E % + -i- + 3^ (^4— J - 2^ ( x + ^ ) + 7s + log(4) ) , (51) 



2V 2 x + 1 T \ 2 J T \ 2 

with H x the Harmonic number, ip(x) the Digamma function and 7^ the Euler's constant 
(7E ~ 0.577 . . .). For j/H ^> 1, the function i? behaves as 



Note we could have derived Eqs. (H7|) and fj49l) working directly in terms of £ = —Hit. 
The only difference would be the replacement N c — >■ iV c = N c /H, and ejv c + e — ► . 
One can also show the terms proportional to 7r not displayed in Eq. (jBJ) , suppressed 
by generalized slow roll parameters, group themselves in such a way to transform the 
equation into an expression as in OH]), but in terms of £ — —Hire with iV c — > N c . We 
already studied some of these terms in the previous section, for more details see appendix lAl 

In the ensuing sections we will demonstrate the validity of Eq. (JTJ at leading non-trivial 
order in the slow roll parameters and mixing with gravity scales. 

3. THE SQUEEZED LIMIT AT FIRST ORDER IN THE GENERALIZED SLOW 

ROLL PARAMETERS 

3.1. Decoupling limit 

In order to set up the stage for the full computation including the mixing with gravity, 
here we start with the ingredients that enter in the consistency condition before adding 
corrections at 0(e), which will be the subject of the next section. In this limit we can work 
directly in the 7r-gauge, since ( ~ — Hit with a constant H. 

Even within this approximation the relationship in Eq. ([1]) becomes non-trivial. For 
example, notice in the equation for ir there is a term fir50 (from expanding f(t + 7r) to 
second order), whose contribution to /nl in the squeezed limit is of order f/(fH) ~ e/ 



111 ] . However, one can imagine a (somewhat tuned) situation in which e Vo + 2e/ 1 



in Eq. ( ]49i) . namely f 2 vo approximately time independent (to preserve scale invariance 
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of the power spectrum) but with a larger Sf, in which case the consistency condition 
might not apply. As we show next that is not the case, and the resolution requires includ- 
ing contributions from the time dependence of the two-point function of the noise, i.e. vo(t). 



We are interested in the three-point function of ( in the squeezed limit. Hence, in addition 
to the non-linear term described above, and the non-linear interaction between n and 50^\ 
we also need the response part SO s R to second order in 7r. Since we are interested in computing 
the non-linear contributions to the equation for 7r in the decoupling limit, namely ignoring 
terms proportional to 7r 2 (and (<97r) 2 ) which do not play a role in the squeezed limit when 
e — > 0, the only piece that contributes from Eq. (fl8l) comes from expanding Vo(t + ir), and 
takes the form 

(50 S R ) {2) (n) - V (t)nn + ■ ■ ■ . (53) 
Collecting all the ingredients the equation for 7r to second order thus reads 

n 2 + (3H + 1 )n 2 -c 2 s ^ = -LtO s + (e Nc - e f ) ^H^SO* + 2s<*Hn 1 ^ 

gt JM C JM C cr 

-e 7 F77Ti7Ti, (54) 



where we only kept terms that contribute in the squeezed limit, and used Eq. 



The source term in Eq. (I54p involving tt50s was analyzed in 11] . The details for the 



computations of /^l for the remaining terms are given in appendix |Bl and we find 

/nl - ~ (2c/ + R(^)e, -3s- 2e Nc ) , (55) 

with R defined in Eq. (15"IT) . Therefore, even ignoring the mixing with gravity, up to this 
point the above expression shows that f^ L is proportional to the tilt if and only if vq is 
constant, namely e„ = in Eq. (|49|) . However, as we anticipated, the key observation is 
that this computation is incomplete, because at non-linear level there is also a contribution 
that arise from the fact that fluctuations in the equal time surfaces, described by 7r, affect 
the probability density functional for the noise when vq depends on time. 

To see how this contribution arises let us exploit the fact that the two-point function of 
the noise is a bi-scalar, hence we need to replace vo(f) vo(t + 7Ti) in the unperturbed 
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expression f[T2"j) 8 

(SOs(x 1 )SO s (x 2 )) = u a (t + ^) 5 ^ Xl ~ X2) ~ */<,(*) (1 + e„ #7n) ^1^1 . (56) 

V 9 a 

This equation suggests (as shown in appendix [C]) one can simply replace: 

50 s <^ 1} (l + ^e^ffTrjJ (57) 

in the equation for tt, while keeping the unperturbed probability density functional for the 
noise. Therefore, the extra piece becomes equivalent to the term which is proportional to Sf, 
up to a factor of 1/2. The upshot of this effect then boils down to replacing e/ — > |e„ + e/ 
in Eq. ( )55l) . This concludes the proof of the consistency condition in the decoupling limit, 
we study the mixing with gravity next. 



3.2. Including the mixing with gravity 

For computing the remaining contributions to /^l that are suppressed by 0(e) we need to 
include the mixing with gravity. Here we will remain in the 7r-gauge. While working with ( 
makes the physics more transparent, the main reason to use tt is that it allows us to proceed 
in complete generality without any particular knowledge about the action for the ADOF. 
That is not the case in the £-gauge, since So depends explicitly on the metric. The price 
to pay is perhaps a somewhat less evident conversion to physical quantities. This, however, 
does not represent a major issue while working perturbatively in the slow roll approximation. 

First of all, notice that terms that involve short-wavelength modes 5N$ and N$ do not 
contribute in the squeezed limit. For starters it is easy to show that terms containing 5Ntt, 
namely without a derivative acting on tt, are multiplied by a background function which is 
already of order one in the generalized slow-roll parameters. This is because factors of tt ap- 
pear after expanding slowly varying functions of time. Since SN is also order one in the slow 
roll parameters, or mixing scales, we can neglect all such terms (see Eq. ( 132]) ). Furthermore 
we can also ignore quadratic terms containing SNstTl- Even though they are not multiplied 



See appendix [C] for a somewhat more formal discussion on the noise kernel. 
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by a generalized slow-roll parameter, in this case ttl yields an additional suppression (recall 
til — cHttl). However, terms involving SNl^s do produce a non- vanishing contribution 
to /nl at leading order. Hence it is only necessary to account for 5N for long-wavelength 



modes. As we showed in sec. 12.31 for the latter we have 5Nl ~ tHu^. 

Taking all these elements into account, expanding the action to cubic order and following 
the standard procedure, the equation for tt to second order becomes: 

7T 2 + (7 + 3H)tt 2 - c\-^- = - (e Nc - 2e) H^fa + 3F7r x ) (58) 

+(2, + e Nc )c 2 s H^^ - L (80f + 50 s ) - jgfHvx (SO® + SO, 

4M 4 Att! 6M 2 4 , 9N 12M 2 4 , 

+ ^T^ - eH ^)^T + -^3^(2eff7ri7ri - tt?) + -j^ir^eH^ - tt x ) 

+ terms that do not contribute in the squeezed limit. 

Recall that the relation between our tt and £ is the same as in single field inflation (which 
is given in Eq. (jUJ)). However, as in the usual case, while in the long- wavelength limit a ( 
mode remains constant, a tt mode does not. Hence a source term tt\ will contribute to /^l- 
It is easy to see using Eq. (jSJ), that 

tt? ~ H + 2eCi% + 0(e 2 ) ~ 2e#7Ti7r 1 + 0(e 2 ), (59) 

and therefore, since ttl — cHttl, the source terms in the last line of Eq. (I5"8~j) do not 
contribute to /^l a ^ leading order. 



In Eq. fl58|) we have both the response and noise parts. For the noise the two-point 
function is given by 

(80 s (x 1 )50s(x 2 )) = u (t + ttO ^ 1 " ^ * »/ (t) (1 + (e„ - e)^) ^ 3 ~ ^ , (60) 

a 6 {l + oN) a 6 

where we used Eq. ( )39|) . and the fact that the two-point function for the noise is a bi-scalar. 
Whereas for the response we ought to use Eq. ( fl8l) . including now the metric components. 
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Expanding to the desired order we get 

{m ,f^_^. m ^.^ ) + ... (61) 

{50 s R f ] ~ d t tttt + S (V - 47V5iV - + (4tt 2 - S^iV) + • ■ ■ 

/jV c7 \ . iV c7 . 
~ a+ — — 7T7T —eHnir + • • ■ , 62 

V / J f 

where in the last line we used eHini + ■ • • , and we only kept terms that contribute 

in the squeezed limit. It is worthwhile emphasizing that, even though naively it seems to 
enter in the computation, the term proportional to Woif) does not enter in the squeezed 
limit after we replace tt — > ns + and SNl — tHir^, to leading order. This represents yet 
another check of the consistency condition, because had Wait) contributed we would wind 
up with a term in the three-point function which is absent in the tilt of the power spectrum. 



Adding all the pieces together the equation for n then becomes 

cr N c JM C a z 

-eHm ( -77T! + 2c 2 ^i - 2^50 s ) , (63) 



c 



where we used the linearized equation for 7Ti, and wrote explicitly the response part SO s R up 
to second order as shown above. 

Using Eq. (|63|) we can compute the different contributions to the squeezed limit of the 
three-point function following the standard procedure. This is done in appendix [B] Finally 
we have to take into account the non-linear relation between ( and tt. Up to 0(e) it takes 
the form ( ~ — Hit + eif 2 7r 2 /2, which corrects /^l by a factor 

/nl I (64) 
Then, garnering all the contributions we obtain (see appendix [B| 

/n! * ~ Uo + 24 - Zs* - 2e% + R (e* + <$) - 2e*) = - 1). (65) 
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One may be worried about the corresponding SA S S , 5B$ from the lapse and shift in 
Eq. (I32|33p . could in principle modify the computation of the three-point function in the 
squeezed limit. That is, however, not the case provided their two-point functions satisfy 
our approximations. The reason is twofold. First of all, they already enter suppressed by 
a power of 1/Mp in Eqs. (132(1331) . Moreover, since they also appear in the same way 5N 
does, terms like 5B s s i[ will be equivalent to those of the form SNir, which as we argued 
before are already multiplied by a background function which is order one, making it higher 
order. Furthermore, since ~ eHni, the only way they could contribute in the squeezed 
limit is via a coupling to "kg. On the other hand, terms like 5B s 7rg can only contribute in 
the squeezed limit from the ir^ dependence in the noise two-point function, for example 
stemming from vs{t + ttl)- However this will not contribute to the three-point function. 
(Also these are suppressed by slow roll parameters, as in Eq. (I6T)|) . thus making their 
contribution a second order effect.) 

This completes the proof of the consistency relation to leading non-trivial order in the 
generalized slow roll parameters. Furthermore, one can also show that the first correction 
to Eq. ( 1651) scales as (ki/ks) 2 , for ki — > 0. This is discussed in detail in appendix IB1 
for the terms we studied in this paper, which contribute to the exact squeezed limit. For 
the remaining terms, such as those involving spatial derivatives of tt, one can show using 
arguments similar to those displayed in j3] that this scaling is not modified. (See appendix 

E) 



4. RESCALING AWAY THE LONG- WAVELENGTH MODE 

In the standard interpretation of the consistency relation 5|-[8j], a crucial point is the 
effect of long-wavelength modes (l on the 2-point function for short scales at zeroth order 
in the gradients of (l, which reduces to a local rescaling of the coordinates £p hys = e^ L x z . 
In this section we show the same occurs in the presence of ADOF, however, for convenience 
we use it rather than (. 9 A long-wavelength mode 7r L considered as part of the background 
corresponds to a deformation of the constant time surfaces of the form t — > t — ttl, so that 



9 Even though, strictly speaking, ttl is not conserved outside the horizon, our procedure is allowed to first 
order in the slow roll parameters, since ~ eHn^. 
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now the perturbations are given by tt = n — = ng- 10 Then the metric in these coordinates 
takes the form 

ds 2 = -dt 2 (l + 5N) 2 + a 2 (t)5 ij (dx i + N i dt)(dx j + N j dt) 

= -dt 2 (l + 5N) 2 + a 2 (t)5 ij e- 2HwL (dx i + N i di)(dx j + N^di) 

= -dP(l + 5N) 2 + a 2 (i)5 ij (dx i + iWt)(dx j + N 3 dt), (66) 

where we defined x l = e~ H7TL x l . Also, to linear order in the perturbations, SN = 5N — -k^ 
and N l = N l (l — n L ). Then, since n L ~ eHn L and 5N L ~ eHn L , we have 

6N L ~ 0, (67) 

hence, to this order, the background metric written in terms of (t, x l ) takes the same form 
as in the original 7r-gauge. Therefore, the equation for its in the presence of hl is nothing 
but a re-scaling of the coordinates: i = t + and x % = e~ HnL x l (or k % = e H7TL k l ). n 
Expanding Eq. ()44l) to linear order in ttl we obtain: 

ti s + (7 + 3#)7r s - c 2 s ^ = - j^SOs + 2eHir L (ns + 3Htt s ) + 2sc 2 s Hir L ^f (68) 

/ 

+jj-(tN c ~ t f )H-K L 50 s + (e - e^Hn L n s 
= ~y 6 °s + (e Nc ~ e f ) -^HttlSOs + 2sc 2 s Htt l ^- 
-e^HTTL-ks + eH7r L I -77i- 5 + 2c^^ - 2^SO s 

where we used 7Tl — ^Ht^l- (To arrive at the second equality we used the linearized equation 
of its-) This is precisely the expression in Eq. (|63|) . 



To conclude our proof an important ingredient is the two-point function for the noise, 



Of course, by construction, the action is invariant under this transformation (provided we also transform 
the metric accordingly). 

Note that the transformation t — > t — ttl, x % — > e~ HnL x l corresponds to an isomctry of the background 
metric at zeroth order in the slow roll parameters. 
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which in the hatted coordinates is given by a similar expression to Eq. ( 1121) . and therefore 
it transforms into: 

(50s(t,k)50s(t'A)) - ^^-( 27 r)^ 3 )(q + k) 

a 3 (t)(l + SN) 

u (t)(l + (e Uo - e)H7r L ) 6 -^^(2n) 3 S^(q + k), (69) 



up to leading order in 7Tl, and we used SNl — 0. This again reproduces Eq. fl60|) . 

5. OTHER TYPE OF COUPLINGS 

In this paper we studied in detail the coupling f(t)0 for the interaction between ir and 
the ADOF. However, there are other possible terms one can write down. We review in what 
follows how to generalize our analysis to those cases. 



5.1. Og 



oo 



The coupling — \Og m was analyzed in detail in ll|. At linear order, and assuming the 
time variation of the ADOF is much faster than the Hubble expansion, namely HSO <C SO, 12 
the only change with respect to the calculation in sec. 12.41 is the replacement SO — > SO, 
together with the requirement ImG® ~ or SO^ ~ 7iV c 7T as in Eq. (S3]), which 

guarantees the 777 dissipative term in the equation for it. Then the computation of the 
power spectrum and tilt follow the same steps. Note however that, by construction, in this 
case the time dependent function / is not present, and therefore one should set ej — > and 
/ — > 1 in the previous calculations. 

Under the local approximation, it is relatively straightforward to show that the results in 
sees. 12.31 and 13.21 still apply, modulo the above replacements. The manipulations, however, 
are somewhat different since the operator —^Og 00 contributes differently to the lapse 
and shift. Nevertheless, one can see that at leading order in the slow roll parameters the 
background part (O) can be absorbed into a renormalization of the p + p term in Eq. ([71). 



Therefore, under the same assumptions as in sec. 12.31 one can show SN ~ eHrr and the 



12 These terms arise after integrating by parts the coupling J d 4 Xy/—ga 3 Oir. 
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constancy of ( holds. 



For the direct computation of the three-point function the non-linear terms in the equation 
for 7r are also not the same as before with a f(t)0 coupling. Moreover, because of the 
requirement ImG^ ~ 1/cu [ll], it is impossible to avoid a non-local behavior in some of the 
non-linear terms that appear in the equation for tt, since under such assumption 50 does not 
have a fully local response to 7r. Similar terms also appear at linear (and non-linear) order 
in the response had we kept the sub-leading pieces proportional to H50 in the equations. 
A proof that includes all such effects is beyond the scope of the present work, however, 
one can show that in the approximation where all such terms are neglected, that is under 
the assumption of a dominant local dynamics, the equation for 7r takes a similar form to 
Eq. (16 3 p (with e/ = 0). The main observation is the following. After absorbing (O) into a 
redefinition of p + p, at non-linear order the contribution from 60 to the equation for 7r can 
be written as 



«r V-9 



-9 pu d p {t + n)d v {t + tt) n^SO, (70) 

where n M is defined in Eq. ( TTTj) . and in the second line we kept only derivatives of SO, as 
explained above. Using this expression, and expanding up to second order we obtain on the 
RHS of the equation for it the terms: 

1 / ~\W 1 / ~\( 2 ) 

-— (l + 7r-SN)(n^O) - — (n^O) (71) 

where we wrote only the terms that may contribute in the squeezed limit and made explicit 
the appearance of the scalar operator: O = n^d^O. Then, we notice that the only difference 
with respect to our analysis in sec. 13.21 is the term: (tc — SN) ( n^d^O J . On the other 
hand, recall tcl — 8Nl ~ 0, therefore this piece does not contribute to the squeezed limit at 
the order we work here when we consider the stochastic part, namely yn^d^dOj , given 
by 80s- However, it may still contribute when instead we take the long- wavelength limit 
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for the response part of the operator and evaluate it on the short mode. But, as in Eqs. 
fl37H6T|) . we have 

(V<V9) (1) ~jN e (ir-5N), (72) 

which as we argued does not enter in f^ L when evaluated for long-wavelength modes. (Note, 
nonetheless, this term produces the dissipative correction for short modes.) This con- 
cludes the proof of the consistency condition for this type of coupling, under the assumption 
of a dominant local interaction. 



5.2. Vectors & Tensors 



For vector and tensor operators the story changes slightly. For example, we can have 
a coupling O^g^ , which may also produce a 777 dissipative term. The difference with the 
Og 00 term is that the former introduces only a linear coupling between d^ir and at the 
level of the action. Again, under the same assumptions as in the previous subsection, the 
computation of the power spectrum and proof of the consistency condition follow similar 
steps as in sees. 12.41 and 13.21 after the identification: 

O -> J=^(V=^n (73) 

Similar considerations apply for tensor couplings. 



6. CONCLUSIONS 



Up until now, checks of the consistency condition had been limited to single field inflation. 
As it is well known, multi-field models with man y lig ht degrees of freedom relevant during 
inflation are capable of violating such condition [24]. In this paper we filled a gap in the 
literature and showed the validity of the consistency relation for a vast class of (multi-field) 
models having a preferred clock. More specifically, using the EFT framework developed in 



Uj , we explicitly demonstrated that: 



lim o (( kl ( k2 ( k3 ) = -(2vr) 3 5 3 (l>) P d k L)Pd k s)(n s - 1), (74) 
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holds for dissipative single-clock inflation, to first non-trivial order in the generalized slow 
roll parameters and mixing with gravity scales. Moreover, we also showed that the first 
correction to the exact squeezed limit scales as (k^/kg) 2 when — > 0. Unlike cases where 
the curvature perturbations are produced in the Bunch-Davies vacuum, here the proper 
treatment of the contribution from the noise to the power spectrum and non-Gaussianities 
played a key role. 



Following our main assumptions were, in addition to the existence of a preferred 
clock, the validity of the local approximation, derivative expansion, and the emergence 
of a shift symmetry at the level of the response for the ADOF. From here we were able 
to show, after including the mixing with gravity, that the effects from the ADOF turn 
off while ( remains constant outside the horizon. The latter was crucial, in particular 
for long- wavelength modes, in allowing us to prove SNl — cHitl in the 7r-gauge, which 



resembles what occurs in single field models. Given the generality of our EFT approach 11] 
we relied on the (non-linear) realization of the symmetries, which we exploited extensively 
in order to incorporate the couplings to tt, as well as the mixing with gravity. Another 
advantage of the EFT formalism, and working with tt even in the long-wavelength limit 
(where the meaning of tt is somewhat less transparent), is that it allows us to bypass any 
specific knowledge about the dynamics of O, other than the analytic properties of Green's 
function and stochastic noise. This would not be the case in the unitary gauge, since ( 
couples to the ADOF in a completely unknown manner in So- 

The results in this paper apply to a plethora of possible scenarios with ADOF, including 
warm [q| and trapped inflation 13 [loj |. and provide further support to the claim that the 
squeezed limit of the three-point function offers a remarkable opportunity to probe the 
very mechanism behind primordial density fluctuations. 14 The forthcoming results from 

as well as large scale structure measurements |l|, |2| and CMB 



the PLANCK satellite 



/x-distortion [3(, thus have the ability to significantly constrain, and/or rule out, a large(r) 
class of inflationary models. 



13 Technically speaking trapped inflation corresponds to a series of operators of the form J2i f&i- It 1S easy 
to show the results of this paper apply to this case as well. 



14 Other soft limits on n-point functions have been recently studied in 26-2 
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Appendix A: C,l in the slow roll approximation 

We want to show that the equation for ( resembles Eq. (I44p to first order in the slow roll 
approximation, including the mixing with gravity. Let us ignore the noise part for simplicity. 
Let us also re- write the background equations in ffl9l - fr2Tl) one more time: 

3H 2 M*=p + p + f(t)0, (Al) 
p+3H(p + p) + f6 = 0, (A2) 
Pa + 3H{p + po) + fO = 0. (A3) 

Expanding the full action the to quadratic order, the linearized equation for 7r (to all 
order in slow-roll) is given by 



1 /rWV 1 ' 

- = -N c {jf - 5N) - {N c + 3HN c )(n - SN) -\p + 3H(p + p)]ir + 3H(p + p)5N 
- f(0 -O)- f~0 • .... (A4) 



(i) 



Na 3 \57T J 



where the dots stand for terms that contain spatial derivatives and N l which we ignore 
since we are interested in the long-wavelength limit where ki — > 0. (Since, as discussed in 



sec. 12.31 the ADOF do not affect the lapse and shift it is easy to show that the terms that 



survive in t 
factors of e 



re long-wavelength limit are proportional to (, and moreover suppressed by 

3 Bo 



Let us start at first order in the slow roll parameters. Note that the equation for it has 
a term which does not involve derivatives, this is required for ( = —Htc to be massless. 
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Using the derivative of Eq. (jA2j) . 5Nl — cHir^ together with ttl ~ eHrci, we find (ignoring 
gradients): 

_L_ (^) {) = £ + 3H(-J^(O-0-dn) = ( + 3H(-i- (50 S R ) {1) (tt), (A5) 

up to second order in slow roll, with N c = N c /H. Hence ( will be conserved outside the 
horizon because 

t (60 R ) (1) (tt l ) ~ jH (TTi - eHn L ) = - 7 Cl, (A6) 

for /cl — ?• 0, which is the case as we argued in sec. 12.31 

We can also demonstrate that the above reasoning applies also to second order in slow 
roll. The main difference, once again, is the appearance of (SOfj)^ = O — O — On — On 2 /2, 
which is then required to be proportional to derivatives of ( in the long- wavelength limit. 
Moreover, all other terms pair up such that the resulting equation for does not have any 
contribution which does not involve derivatives. 



Appendix B: The source terms in the squeezed limit &; {ki/ks) 2 scaling 

In this appendix we provide some details on the computation of /^l f° r the different 
types of source terms appearing in Eqs. ( J54l) and (l63~l) . We also show that the first 
non-trivial correction in the squeezed limit scales as (ki/ks) 2 - 

Let us start with the source term that is proportional to €f\ —N~ x ftfHndOs- The 
contributions of the other terms involving nSOg can be obtained from here by performing 
the appropriate replacements. To analyze the non-Gaussianities we decompose n = i\\ + n 2 , 
where again the subscripts represent the order of the solution for the fluctuations. Then 

fkc s P 

tti(M) = j d v'9i(kc a \7i\,kc s \rf\)50s, (Bl) 

where g^(kc s \rj\, kc s \rj'\) = G~,(kc s \ri\, kc s \rj'\)/ (kc s rj') 2 with G 1 defined in Eq. (HEJ), and 
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(using Eq. ([EI]) ) 

, e f f 2 k 3 c 2 f° . . f rf 3 q 

Ws,0) = N2H3 J drjg^ksCslri |) y 

x f rfr / "( ?7 (g Cs |r / '|,g Cs |r / "|)505(q,V / )^5(k 3 -q,r/ / ). (B2) 
^ no 

We want to compute the three-point function for ( ~ — i?7T in the limit i] — > 0: 

(C(ki, 0)C(k 2 , 0)C(k 3 , 0)) = (d(ki, 0)Ci(k 2 , 0)C 2 (k 3 , 0)) + cyclic sum in fc's 

= (2ir) a S 3 (^2k i )F{k 1 ,k 2 ,k 3 ). (B3) 

i 

Then, using Eq. ( |T2l) and that at leading order the noise is Gaussian we have 



(SOsik^^SOs^fj^Osiq^yOsiks - q, r/")> = (2vr) 6 4 Q 4 ( ^ >4 ^ w) 7 (B4) 
{<%' - rf')6{f} - 7/")^(k 2 + q) + 5{fi - r}")5(fi - 7/")<* (3) (ki + q)} (for A; 3 ^ 0). 



x 



Defining Xi = fcj/A;, with /c an arbitrary scale with units of momentum, we obtain (for 
i] -oo): 

F(xi, x 2 , x 3 ) = k 6 F(h, k 2 , k 3 ) = ° J \ x x\x$ \ dy y^g 7 (0, x 3 y)gj(0, x x y) 

p+oa 

x / dz z 4 g~ / (x 2 y, x 2 z)g 1 (0, x 2 z) + permutations in x^s (B5) 
Jy 

(we performed a change of variables: y' = —kc s T]', z = —kc s rj"). 

From the definition of /^l gi ven i n Eq. (jlj) and the expression in Eq. f )47|) for the power 
spectrum, it is then straightforward to show 



/n1| ="gC/, (B6) 



independently of the value of 7. 



We follow now a similar procedure for the contribution due to the change in the probability 
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density functional of the noise. In this case we need 



(C(k 1 ,o)c(k 2 ,o)c(k 3 ,o)) = - /3 ^!f 3 f dm f d V2 f ^ 7 (o,cA|^i|)(B7) 

■'"c J —oo J —oo J —oo 

x# 7 (0, c s k 2 \r] 2 \)g^(0, c s k 3 \r) 3 \){SO s {m, ^i)SO s (v2, k 2 )50 s (r] 3 , k 3 )). 



Using Eq. (JDZD and 



(SOs^m^OsQLi, V2 )SO s (k 3 , V3 )50 s (k, v')SO s (p, r]')80 s {-k - p, r/')) = (B8) 

= (2 " )9 ^ S^m^ - ^ - ^ - ^ ,)5(3)(ki + k)5(3)(ks + p) 

+ permutations in the fcj's} , 

one can show that this term yields a contribution to F(xi,X2,x 3 ) that is exactly the same 
as the one in Eq. (1B5j) after replacing e/ by e vo /2. 



Let us now analyze the contribution to /^l from the last term on the RHS of Eq. (l54l) : 
—e^H'-fTTTc. This term contains a time derivative of it, but it also involves an additional 
temporal integration. For this reason it is non-trivial to see whether this vanishes or not in 
the squeezed limit. In this case we have 

e if 3 7^ 2 f A [ + °° f +co f +co 

F(x 1 ,x 2 ,x 3 ) = 7 6 ° — x x x\x\ / dy dz dwyz A w A g 1 (Q ) xiy) (B9) 

^c c s JO Jy Jy 

<? 7 (0, x 3 z)g~ f (x 3 y, x 3 z)g y (0, z 2 w)— g^(x 2 y, x 2 w) + permutations in x^s. 

ay 

In the squeezed limit x 3 — > 0, X\ — > x 2 = 1, the contribution that does not vanish is given 

by 

2e if 3 7z/ 2 / 4 /" + °° /" + °° /" + °° 
C= Kma;lF(l,l,Z3) = J™ 7 n a c 6 J q d V J d ~ z dwyz 4 w%(0,y) 

flf 7 (0, z)g^(x 3 y, % 7 (0, w)-^g^{y, w) 



2e if 3 7z/ 2 / 4 /" +oc /" +oc /" + °° 

L ' / c?y / / dwyz^w^g^O^y) 



flf 7 (0, % 7 (0, % 7 (0, w^g^y, w) (BIO) 
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where we made the change of variables y = x±y, z = x 3 z, w = x 2 w and we set xi — x 2 — 1. 
The contribution to /^l rea ds 

'»\=^i#r-^ B (tf)' (B11) 

with R(y/H) given in Eq. ( |5T|) . 15 

Similarly we obtain a piece proportional to e from Eq. ( 1631) after replacing e 7 — > e, and 
for the source term 2sc s Hir 1 ^ L we find: /^i = |s, independently of 7. 

I s 

To probe the scaling with ki/ks from each one of these terms it is sufficient to take 
each contribution and expand in powers of x 3 . For suggestive purposes in what follows 
we re-label 23 — > xl- To show that the first non-trivial correction starts at 0(x 2 L ) we just 
need to show that the first derivative of x 3 l F(xs,xs,xl) vanishes at xl = 0. Here we 
assume analyticity of the result in xl- (This is not required by the consistency condition, 
since for instance we could still have y/xZ- However, by inspection it is easy to show each 
contribution is indeed analytic.) 



As an example, let us consider the term in Eq. (IBlOp but for non-zero xl- If we take 
now dC g*^ 1 we encounter a few different contributions. First of all we have the piece that 
depends on the limit of integration, and therefore we need to evaluate the integrand at xl- 
This produces pieces proportional to g 7 (0, 0) which vanish. The remaining term depends on 



the integral o: 
the property 



a derivative of the Green's function, for example d XL g-y(xiy, z). Then, using 
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d XL g 7 (x L y, z) ~ — — —y g^(x L y, z), (B12) 
(1 + J/H) 



which is valid for yx^ <C y* — yy/H, and taking into account that the contribution from 
g 7 (0,?/) in the integral is dominated for values y ~ y*, 16 one concludes that this term scales 
like Xl for xl 1, hence J^- vanishes when xl — > as advertised. 



15 Even though this equality is not straightforward, one can show that the integrand in Eq. (|B10| agrees 
with an equivalent representation of Eq. (|5ip term by term in an analytic expansion, once the common 
non-analytic pieces^yie factored out. Moreover, the results also agree numerically. 

16 See appendix F in 11] 
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One may wonder whether the terms we ignored in this paper, that do not contribute 
in the exact squeezed limit, could in principle produce a correction at order k^jks- As 



shown in that is not the case for the pieces proportional to 50ir, since in fact they are 
sub-dominant in the squeezed limit. Moreover, following one can also show that the 



other possible terms studied in 11], such as ^{dir) 2 etc., do not modify the scaling in Eq. 
([I]). The basic idea is the following. Since the lapse and shift are determined by the same 
expressions as in single field models, the arguments put forward in [7] translate directly to 
our case. For example, one can have a term of the form N l di7r50, from the mixing with 
gravity. But SN ~ and N l ~ O(ki), and therefore there is no left over k^/ks upon 

symmetrization [7|. Note also that terms coming from the vector coupling O^d^Tc turn into 
a scalar interaction of the form (d^O^ir. Moreover, couplings of the sort g^O^ do not 
produce non- linear terms in ir other than through the non-linear response for d^O^, which 
is included in our analysis. 

Appendix C: Non-linear generalization of the probability density functional for the 

noise 

The issue of including the effect of ir on the probability density functional for the noise 
can be naturally addressed in the framework of the closed time path (CTP) or Schwinger- 
Keldysh formalism, where one has a CTP effective action for the "mean" field, that is given 
by 7r + = (7T 1 + 7r 2 )/2, and the difference 7r_ = 7T 1 — 7r 2 (in CTP notation). In general, the 



effective action has real and imaginary parts (see for instance |29|). 



exp{— iT[n + , 7r_]} = exp{— ?3?r[7r + , 7r_] + Sr[7r + , 7r_]}. (CI) 

The imaginary part can be interpreted as the result of an averaging over an stochastic source 
£ = 60s. For example, for the coupling f(t)0 and to second order in tt, the imaginary part 
of the effective action takes the form 

eXp{QT}[7T + , 7T_] = F[7r_(x)] 

= exp j-^y d^a^h) J d 4 x 2 a 3 {t 2 ) A_(x 1 )A^(xi,x 2 )A_(x 2 )| , (C2) 
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with A_ = fn 2 + f(n 2 ) 2 /2 — fn 1 — f(n 1 ) 2 /2, and N(xi,x 2 ) is the noise kernel. The prob- 
ability density functional of the stochastic source is then obtained after using the following 
mathematical identity 



F[ir-(x)] = J 'z>f(x)P[£(a;)]exp j-i J d 4 ya 3 A_ 



-{V)Z{V)}, (C3) 



with 

ft rt 



P[Z]=Aexp[~ J ^jo 3 ^) J' d 4 x 2 a 3 (t 2 ) ^x 1 )N- 1 (x 1 ,x 2 )^(x 2 )Y (C4) 

where A is a normalization constant and N~ 1 (xi,x 2 ) is the inverse functional of N(x\,x 2 ). 
In the local approximation (in our case) 

N(x 1 ,x 2 ) = u {t f {Xl ~ X2) (C5) 

cr 

iV (zi,^) = — 7- ; • (C6) 

v (t) a A 

The equation for n, as an stochastic variable, is then obtained before averaging over £, 
namely from the real part of the action where we have the term J t d A y a 3 A_ (?/)£(?/). 17 The 
contribution for this term is equivalent to our fSOs + fefHndOs in Eq. ( 1541) . 

In cases where the probability distribution depends on n + (i.e., P[£] — >■ P[£, 7r+]) 
the previous identity remains valid. This justifies the replacement in the noise kernel 
vo —> vo{t + 7r) (recall that at the level of the equation of motion n + = ir). 



Our objective is to compute (nun) perturbatively. So, we split n = n\ + n 2 , where n 2 
(like 7Ti) has 50s = £ as a source, but for n 2 the averaging over the noise (. . .) is computed 
using P[£, 7Ti], which is obtained after the replacement vq — > vq + ^o 71 "!- Therefore, we need 
to know obtained using P[£, 7Ti]: 

(eee) * | T>tP[t,*+ = o] Q 1 1 dVa^odW^eiv-^Trx) eee[ =o - (cwo 

Now we have to replace n% by the solution in Eq. ( 1431) . Note that there are six stochastic 



17 The variation can be taken either with respect to ir 1 or ir 2 , and afterwards one sets n 1 — tt 2 (i.e., 7r+ = 7r 
and 7T_ =0). 
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fields (£'s) in Eq. (10T|) but in the denominator there are two factors of uq (one is inside N) 
and we have a factor 1/2. It is not difficult to see that the result is exactly the same as the 
one for —N^fefHn^Os (see appendix IB1) . but replacing ey — » e vo /2. Then, the effect of 
both contributions becomes 

/n!| ^ — (6^+26,). (C8) 

Perhaps a more straightforward way to show that this is the combination that appears 
in the three-point function is to use A_ = fn_ + fn_n + , and rewrite Eq.( 1C3j) as 

F[n-(x)] = J P£(x)P[£(x)]exp{-* j* d 4 ya 3 fir-{v)^v) } , (G9) 

where in the equation for 7r 2 (Eq. ([51]) ) we replace the source term f50$ + fejHiTi50s by 
fSOs, but with a different probability density functional: 

P[£] = Aexp|-iy rf 4 xia 3 (ti) y rf 4 x 2 a 3 (t 2 )e(xi)iV- 1 (xi,x 2 )e(x 2 )| , (CIO) 



with 



a 3 



z^ 1 (1 - {e vo + 2 e/ )tfvr + ) ( C ll) 

CL 



This analysis provides a formal support to the manipulations behind Eqs. (16T)1) and (]69[) . 
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